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Spin dependent electron transport through a magnetic resonant tunnehng diode 
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Electron transport properties in nanostructures can be modeled, for example, by using the semi- 
classical Wigner formalism or the quantum mechanical Green's functions formalism. We compare 
the performance and the results of these methods in the case of magnetic resonant-tunneling diodes. 
We have implemented the two methods within the self-consistent spin-density-functional theory. 
Our numerical implementation of the Wigner formalism is based on the finite-difference scheme 
whereas for the Green's function formalism the finite-element method is used. As a specific appli- 
cation, we consider the device studied by Slobodskyy et al. [Phys. Rev. Lett. 90, 246601 (2003)] 
and analyze their experimental results. The Wigner and Green's functions formalisms give similar 
electron densities and potentials but, surprisingly, the former method requires much more computer 
resources in order to obtain numerically accurate results for currents. Both of the formalisms can 
successfully be used to model magnetic resonant tunneling diode structures. 

PACS numbers: 73.63.-b,75.47.-m 



I. INTRODUCTION 

The future spintronics technology requires controlled 
spin injection into semiconductor materials. The prob- 
lem can be solved using different techniques and con- 
structions (see, for example, Refs. 2, 3, ^, ^^^J). The 
magnetic resonant tunneling diode (RTD) demonstrated 
by Slobodskyy et al. [J is one of the most promising 
solutions. Their magnetic RTD is based on the quantum 
well made of dilute magnetic semiconductor ZnMnSe be- 
tween two ZnBeSe barriers and surrounded by highly n- 
type ZnSe. In the presence of a magnetic field there is 
a giant Zeeman splitting between the spin-up and spin- 
down electron states in the quantum well region. The 
corresponding resonance peaks in the electron current 
separate as a function of the bias voltage, and the de- 
vice controls the electron spin states using the bias volt- 
age rather than external magnetic fields. This is a useful 
property for possible practical applications. 

In a typical nano-scale transport problem two or more 
electrodes are connected to a functional nanostructure. 
An important nano-system is the RTD composed of lay- 
ers of semiconducting materials. The offset between the 
band edges of the materials give rise to the two potential 
barriers seen by carriers. In the quantum well between 
the barriers, resonance states with finite energy width 
and enhanced amplitude are formed. The electron cur- 
rent increases rapidly when a resonance state appears 
in the conducting window, determined by the overlap 
of the occupied source and unoccupied drain conduction 
electron bands. When the resonance state drops with in- 
creasing bias voltage below the source conduction band, 
the current diminishes causing a region of negative resis- 
tance. Besides being technologically interesting the RTD 
is also important as a simple test case for different for- 
malisms and computational schemes. 

Theoretical modeling and computational simulation 
are essential for the development of functional nanos- 
tructures. Electron transport properties can be modeled 



using different formalisms at different levels of sophistica- 
tion. Two methods widely used are the Wigner function 
(WF) 0,13 and the Green's functions (GF) formalism 0]. 
The WF approach is a semiclassical transport formalism 
that enables the study of systems exhibiting quantum in- 
terference and tunneling effects. The electron density and 
the current are obtained from the Wigner function, which 
is in turn calculated by solving the Liouville-von Neu- 
mann equation. The GF formalism is a fully quantum 
mechanical scheme with an increased complexity with 
respect to the WF formalism. Both the WF and GF 
formalisms enable self-consistent non-equilibrium calcu- 
lations corresponding to a finite bias voltage between the 
electrodes. 

RTD's can be modeled as structures which are trans- 
lationally invariant parallel to the layers. This makes the 
computations one-dimensional. The WF and GF for- 
malisms are very popular schemes in their modeling (see, 
for example, 0, 0, 0|). In this work we study 
the feasibility of the WF and GF formalisms within the 
context of the self-consistent spin-density-functional the- 
ory (SDFT) and its local spin-density approximation to 
model magnetic RTD's. Our implementation of the WF 
formalism is based on the usual discretation of the partial 
differential equations on position and momentum point 
grids. Our implementation of the GF formalism employs 
the finite-element method (FEM). We have already pub- 
lished our FEM scheme for two-dimensional nanostruc- 
tures 0] and used it in applications . In this work we 
critically compare the performance and the results of the 
WF and GF implementations, for quasi-one-dimcnsional 
RTDs. Moreover, we choose the structure parameters of 
our test system to correspond the magnetic RTD device 
by Slobodskyy et al. Thus, we can compare our 

results also with experiments and actually analyze the 
results of the measurements. 

Below we use the effective atomic units which are de- 
rived by setting the fundamental constants e — h — rrie — 
1, and the material constants to* = e = 1. m* and e 
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are, respectively, the relative effective electron mass and 
the relative dielectric constant to be used in the effective 
mass approximation. For ZnSe m* = 0.145 and e — 9.1 
[TgL [itII . We have used same values also for the other ma- 
terials of the magnetic RTD. The effective atomic units 
can be transformed to the usual atomic and SI units using 
the relations 

Length: 1 = ~ 3.32 nm 

Energy: 1 Ha* = 1^1^ Ha « 47.6 meV 
Current: la.u.* = l^a.u. « 11.6 fiA. 

Above, Ha denotes the Hartree energy unit. 

In Sec. II we introduce the model for the magnetic 
RTD and explain the use of the SDFT in the calcula- 
tions. In Sec. HI we briefly explain the GF and WF 
formalisms. In Sec. IV we give the results of the com- 
parisons between the two formalisms and compare the 
calculated and measured results and obtain information 
on the electronic structure of the device in question. Sec. 
V contains our conclusions. 



II. MODEL 
A. Structure of magnetic RTD 

The model for the magnetic RTD is shown in Fig. ^ 
We assume the semiconducting layers to be infinitely 
wide in the lateral directions so that the system is trans- 
lationally invariant in the direction perpendicular to the 
current. The doped regions (shadowed areas in Fig. ^ 
are modeled by a uniform positive background charge. 
The potential barriers due to the discontinuity of the 
conduction band between two materials are described by 
constant external potentials. 

The system is divided into the central region f2 and 
the outside regions fi^ and flR. All the structural vari- 
ations and interesting phenomena take place in region f2 
which is chosen large enough so that the effect of the RTD 
device on the electron density has vanished at the bound- 
aries dVLj^/n. ^L/n are the semi- infinite leads where the 
electron density and the potential are constant. The bot- 
tom of the conduction band and the Fermi- level in the 
right lead are shifted by the bias voltage Vsd with respect 
to the corresponding values in the left lead (see Fig. 
In our model the electron transport is ballistic with no 
phonon or defect scattering. This means that the total 
potential drop takes place within Vl. 

In addition, f7 is divided into nine smaller parts SI, ... 
S9 as shown in Fig. ^ These regions describe different 
semiconductor material layers. The parameters of the 
layers are given in Tabled We have chosen the widths 
and the doping densities of our magnetic RTD structure 
similar to those in the actual device made by Slobodskyy 
et al. The positive background charge densities in 
regions SI and S9 are equal to that in fii/^j, in S2 and S8 
they are much smaller, and next to the potential barriers 
there is no background charge at all. 
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FIG. 1: Magnetic RTD model. The shaded areas denote 
positive background charge. The detailed information about 
the different layers S1,...,S9 are given in Table in 



The regions S4 and S6 are the potential barriers. Be- 
cause there is no definitive information about the barrier 
height we estimate it by comparing the shapes of the cal- 
culated current-voltage curves to the measured ones. By 
the shape we mean mainly the widths of the resonances; 
the positions of the resonances are quite insensitive to 
the barrier height. The barrier height is the only struc- 
tural parameter which we have to determine by fitting. 
We find that the barrier height of 23% of the band gap 
difference between ZnSe and ZnBeSe results in a 

good fit. 

The quantum well S5 is made of the dilute mag- 
netic semiconductor ZnMnSe. An external magnetic field 
causes a giant Zeeman effect, splitting the spin- up and 
spin-down electron states. In the other parts of the RTD 
the spin splitting is small and assumed to vanish. We 
calculate the spin splitting Ai? in S5 as Slobodskyy et 
al. Jj, i.e. 



AE = NnaxsoB, 



\kB{T + T,ff) 



(1) 



Above, Nqu is the s-d exchange integral, x the Mn con- 
centration, g is the Lande factor, hb the Bohr magneton, 
Bs is the Brillouin function of spin s, sq is the effective 
Mn spin and Te// is the effective temperature. The val- 
ues of the parameters are N^a — 0.26 eV, x — 8%, s = 



TABLE I: Parameters of the different layers SI, 
the calculations. 



,S9 used in 



Region 


Material 


Width 
(nm) 


Doping level 
(cm-3) 


Potential Vw 
(meV) 


SI, S9 


ZnSe 


25 


15x10'** 





S2, S8 


Zno.97Beo.03Se 


15 


1x10^** 





S3, S7 


ZnSe(i) 


10 








S4, S6 


Zno.7Beo.3Se 


5 





92 


S5 


Zno.geMno.osSe 


9 
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TABLE II: Values of the spin splitting /S.E in the cases con- 
sidered in Sec. V. 



Magnetic field (T) 


Temperature (K) 


AS (meV) 





4.2 


0.0 


2 


4.2 


10.4 


4 


4.2 


16.7 


6 


4.2 


19.8 


6 





23.2 


6 


8 


16.2 


6 


30 


6.7 



5/2, g = 2.00, so = 1-13 and Te// = 2.24 K. The values of 
the AE for the relevant magnetic fields and temperatures 
(Sec. V) are collected into Table Hll 



We calculate the currents due to the spin-up and spin- 
down electrons through the magnetic RTD as a function 
of the voltage. The calculation always starts from the 
zero bias voltage Vsd- When the self-consistent non- 
biased (equilibrium) result is reached, we increase Vsd hi 
small steps and iterate at every value until convergence. 
The effective potential corresponding to the previous Vs d 
value is used as the starting point of the iterations. This 
ensures the stability of the process. Actually we found 
that the resulting current voltage curve does not depend 
on the sweep direction of the bias voltage. This is due to 
the very small electron density in the quantum well region 
S5 even in the case of occupied resonance states. This 
insensitivity is in accord with the findings by Slobodskyy 
et al. 



B. Spin-density- functional theory 

In order to model electron-electron interactions we use 
the SDFT within the local-density approximation. The 
electronic structures and currents are calculated using 
the GF or the WF formalisms, which we will explain later 
in Sec. IIL In both formalisms the spin-up (cr =t) and 
spin-down (cr =J,) electron densities n-f^|(x) correspond 
to the effective potentials 

K'} f (x) = Vc {x) + V^, (x) + (x) , (2) 

where Vc is the Coulomb potential, V^^ the exchange- 
correlation potential |^ and V^ the external potential 
including the barriers and the giant Zeeman spin splitting 
/\E (see Table inil. The spin densities and the effective 
potentials are solved self-consistently. 

The Coulomb potential Vc is calculated using the mod- 
ified Poisson equation [23| 

V^V^+\x)~eV^+\x) 

= - 47r[n+(x) - n'_{x)] - k^V^{x), 

where n_(x) = n-^{x) + ni{x) is the total electron den- 
sity and n+(x) is the positive background charge. Index 
i counts the self-consistency iterations, so that V^{x) is 
the solution from the previous self-consistency iteration. 
Above, fc is a parameter which controls the screening of 
the potential fiuctuations due to the charge sloshing be- 
tween the iterations. A reasonable choice of the k param- 
eter is of the order of the Thomas-Fermi wave vector, in 
which case the solution does not depend on k and the 
number of the self-consistency iterations needed is re- 
markably reduced. Besides, the stability obtained by the 
use of the modified Poisson equation, we stabilize the it- 
erations also by mixing the old effective potential V^jj: 
with that obtained from Eqs. (|2Jl and J^Jl. I.e., 

V:+' =aVff + il-a)V:jf, (4) 

where the feedback parameter a is typically 0.2 in our 
calculations. 



III. FORMALISMS 

In this section we present the GF and WF formalisms 
by ignoring the spin-dependence for simplicity. The gen- 
eralization to the spin-dependent forms, which we use in 
the actual calculations, is straightforward. 

A. Green's function formalism 

The GF formalism used in the electron density and 
transport calculations is explained in detail in Ref. |^. 
We have implemented this formalism using the finite el- 
ement method (FEM). Our FEM formulation for two- 
dimensional nanostructures is discussed in Ref. |l4l| . 

Our present magnetic RTD system is translationally 
invariant in directions perpendicular to the electron cur- 
rent. The total electron energy ujtot can then be divided 
into two parts 

wtot =uj + uj±, (5) 

where lo± is the kinetic energy in the perpendicular direc- 
tions and Lo includes the kinetic energy along the current 
and the (one-dimensional) potential energy. Now we can 
write a one-dimensional equation for the single-particle 
retarded Green's function in the real space using spatial 
coordinates, x and x' in the direction perpendicular to 
the layers: 

{uj~H)G''{x,x';uj)=S{x~x'), (6) 
where H is the Hamiltonian 

H = -l\7' + V,ff{x). (7) 

Above, uj has a small imaginary part, i.e., to — lo' + irj. 
This distinguishes between the retarded and the ad- 
vanced Green's functions. Eq. © is solved using open 
boundary conditions at dn^ and dfln. This means that 
electron wave functions penetrate the boundaries without 
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reflection. In order to restrict the numerical calculations 
into the central region Q Eq. (jSJ is written in the form 

{u;~Ho-^1{u;)~J:''j,{u;))G^{x,x';lj)^5{x~x'), (8) 

where ^-re the self-energies of the leads ^l/r 

and Hq is the Hamiltonian of the isolated region fl. 

The electron density is obtained from the so-called 
lesser Green function by integrating over uj 



-1 f°° 

n-{x) = — / lra{G^ {x, x;Lu))duj. 
27r J_oo 



(9) 



When no bias voltage is applied the system is in equilib- 
rium and is calculated as 



G<{x,x';u;) = 2Fl/r{w)G^x,x' -uo) 



(10) 



Above, -Fl/jj are related to the Fermi distributions in 
^L/B. Ei-nd in the equilibrium they are equal. The effects 
of the perpendicular directions in our computationally 
one-dimensional system are included in Fj^/n by inte- 
grating the Fermi distributions over LOi_ 



-ksTln I 1 + e 

TT 



■duj± 



(11) 



When the bias voltage Vsd is appHed Fl{uj) and Fji{uj) 
are split by Vsd on the energy axis. In this case has 
to be calculated as 

G^(x, x'; uj) = 

- iFR{uj)G'^{x, xr\ uj)Tr{xr, x'r\ uj)G'^{x'ji, x'; uj) 

- iFL{uj)G''{x,XL;LL>)TLixL,x'^;uj)G"'{x'^,x';uj), 

(12) 

where xi^ /r are the coordinates of the boundaries 517^ /r 
and T^/R are defined as 



iT 



L/R - ^L/R - ^L/R - 2«Im(Si/^). 



(13) 



Eq. (|12|l is valid also in equilibrium, but only when there 
are no bound states. Namely, in this form the electron 
density in f2 is composed of scattering electron states 
coming from f2i and Qr. 

To perform the energy integral in Eq. is the heav- 
iest part of the calculations. The calculation of G^ at 
several energies lo takes CPU time but cannot be avoided 
due to the sharp resonance peaks in our RTD system. 
In order to reduce the number of uj values needed we 
move parts of the integral to the complex plane where 
the changes in G^ are smooth. The move of the inte- 
gration part away from the real axis requires that the 
integrand is analytic above the real axis. To fulfill this 
we first approximate Eq. Hll(l as 




Re(co) 



Re(co) 



FIG. 2: Integration paths for the electron density calcu- 
lations. The paths are divided into three parts n_i, n_2 
and a) and b) correspond to the cases Vsd < Hl and 

Vsd > HL, respectively. 



This is exactly true at zero temperature and a good ap- 
proximation at energies few fc^T below the Fermi-levels 
Mi/fl- Now we can split the integral into three parts n-i, 
n_2, and n_3 shown in Fig. [51 The slightly complicated 
division is due to the proper inclusion of the perpendicu- 
lar kinetic energy component. Here and below we assume 
that fiL < ^J'R■ In the first part n_i the form of G^ in 
Eq. Hl()(l is used and the integral 



-(/Xi — w)Im [G''(x, x'; (15) 



is moved to the complex plane. Above, A is about 3 ... 
5 so that Eq. H14I) is approximately valid. The integral 
starts at the energy wq, which is below the bottom of the 
conduction band of the left lead so that the contribution 
of possible bound states is also included. The integral 
ends at a couple of fc^T below the Fermi-level ^l. 

The integrals n_2 and n_3 do not vanish in non- 
equilibrium. Because Eq. (|12|l has to be used the in- 
tegrands are not analytic on the complex plane and the 
integrations have to be performed along the real axis. 
The integral n_2 does not vanish when Vsd is smaller 
than the width of the occupied conduction band in the 
leads iR (the case of Fig. |2i) and is obtained as 



n^2 = I -(Mi? - AiL)x 



TT 



(16) 



G'-ix, XR)rR{xR, x'R)G''ix'j,, x')dw. 



(14) 



The integral n_3 starts at a couple of /cbT below /i^ and 
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ends at a couple oi ksT above /ijj. It is obtained as 



n-3 = 



1 

2^ 



[ Fr{uj) ^(x, xr) Tr{xr, x'r) G'^ix'R, x') 

+ Fl (c^) C {x,xl) Tl {xl ,x'j^)G'' {x'^ ,x')\dLo. 

(17) 

When VsD is larger than the width of the occupied con- 
duction band in the leads ^hjR (the case of Fig. IJb) 
electron states from /i^ + ^ksT to Vsd are not occu- 
pied. 

The electron tunneling probability through f2 is calcu- 
lated using the Green's functions as 



T{uj) =rL{xL,x'i^;u;)G''{x'^,xii;u;) 



><^R{xR,x'jf-uj) G°-{x'j^,xl;uj) 
and the current is obtained by integrating over w, i.e. 
1 



(18) 



TiLu){flLH- flR{L0))di0. (19) 



This whole integral has to be calculated on the real axis, 
which is not a problem because the integral is evaluated 
only once after the self-consistent solution of the electron 
density is reached. 



B. Finite-element implementation of the GF 
formalism 

We have implemented the GF formalism using the 
FEM. For the FEM implementation we need to write 
the equations to be solved in the variational form. We 
start from Eq. (|SJ| and multiply both sides by a contin- 
uous well-behaving function v{x) and integrate over Q. 
After modifications the equation obtains the form 



-Vv{x) ■ -WG'\x,x';lu) 
+ v{x) [uj — Vey/(x)] G"'(x, x'; (jj)| dx 

+ v{xi}j i:l{xL,XL) G^'ixL.XL-^Uj) 

+ v{xr) t'RixR, xr) G''{xr, xr; Ll>) 
= v{x'). 



(20) 



Here the self-energy operators S^^^ have the analytic 
solutions Q 



I ,,1^''91/r{x.x'-Uj) 



(21) 



dxdx' 

where ^'"^ retarded Green's functions in the iso- 
lated leads ^L/R so that they vanish at the boundaries 



an 



L/R- 



In our model the potential is constant in the 



leads so that 



^iV2Lo(x-x') _ i\/2iZ(x+x' -2x L / a) 



2UJ 



In the FEM the retarded Green's function is expanded 
in the basis (bi. 



N 

G''{x,x') « ^ gij4ij{x)(j)i{x'). 



(23) 



This is an approximation due to the finite number N of 
the basis functions. Eq. (|23|l is inserted into Eq. H20|l . 
The values of the coefficients gij are then calculated by 
choosing v{x) = (l)k{x). 

We use a basis consisting of the linear functions (fP and 
(jx^ and of high-order polynomials ,2^. That is, 



= ^ 2(2/- 1) ^■'■-2 ' ^ 



= 2,3, 



(24) 

where Pj(x)'s are the Legendre functions of the or- 
der j. '/^'(x)'s are given in the reference element with 
X = [—1, 1]. The linear functions span the region of two 
elements, whereas the j > 1 functions are localized within 
one element only. The inclusion of the high-order basis 
functions reduces the number of basis functions needed in 
order to obtain accurate results. This has a remarkable 
effect in GF electron structure calculations even for one- 
dimensional systems. The calculation of the G""(a;, x'; w) 
requires the inversion of a matrix of the size oi N x N . 
The derivatives of the functions with j > 1 are orthog- 
onal to each other. This makes their use numerically 
stable and we have implemented elements up to the fifth 
order. 



C. Wigner function formalism 

We have implemented also the WF formalism for the 
electron density and the current calculations. Reviews 
of the WF formalism for studies of open systems can be 
found in Refs. and There are similarities between 
our GF and WF implementations in the treatment of 
the one-dimensional equations resulting from the trans- 
lational invariance in the directions perpendicular to the 
electron current. They show up, for example, in the 
boundary conditions as will be discussed below. 

The derivation 0, Q of the WF formalism starts 
by making a coordinate transformation to the quantum- 
mechanical density matrix 



p{x,x') = Wi{x\i){i\x') 



(25) 



(22) 



where x and x' are spatial coordinates, |i) is a complete 
set of states, and Wi gives the probability of finding a 
particle in state The new coordinates, q and r, are 
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jiven by the relations 



and 



q = -{x + x) 
r = {x — x') . 



(26) 



The WF f{q,p) is now defined as the Fourier transform 
of the density matrix, i.e. 



fip,q) 



-ipr 



1 



1 



p{q + -r, q ~ -r)dr (27) 



and a classical phase-space representation f{q,p) is ob- 
tained. 

In order to use the WF formalism in transport the- 
ory we must study the time evolution of the WF. Taking 
the time derivative of the density matrix and substitut- 
ing it into the Schrodinger equation gives the quantum- 
mechanical Liouville-von Neumann equation 



■^P !u 1 r 
I'-qI [H-,p\ = Cp, 



(28) 



n-{q) 



27r 



f{q,p)dp. 



(33) 



In an open system we need to make a difference be- 
tween the incoming and outgoing particles to account 
for the irreversibility. Our system is one-dimensional, so 
that < q < Wn, and there are two boundaries at g = 
and q = Wq. Here, Wq = xr — xl '^s the width of the 
central region fl. Since the characteristics are first-order 
equations only one boundary value is needed. Moreover, 
since particles which have p > are moving in the posi- 
tive direction on the q-axis, we must supply the bound- 
ary conditions on the left hand side boundary from 
where they are originating with the momentum distribu- 
tion in the reservoir. Similarly particles with p < are 
moving in the negative direction on the g-axis so that the 
boundary conditions must be specified on the right hand 
side boundary dflR. This scheme invokes the boundary 
conditions 



where H is the Hamiltonian of the system and C is the 
Liouville superoperator. Substituting the Hamiltonian H 
of Eq. into Eq. |2HI) , we obtain 



.dp _ 1 f 



P 



[Veffix) 



dx'"^ 

Veffix')] p, 



(29) 



where Vg// is the SDFT effective potential given in 
Eq. If a Fourier transform, similar to that for the 
density matrix, is performed on Eq. (|29|l the transport 
equation for the WF becomes 



df{q,p) df{q,p) 
= -P- 



dt 



dq 
' l_ 
„ 2^ 



V{q,p~p')f{q,p')dp' .m 



Above, 



f{Q,p%>O^FL ( \p^ 

f{Wn.p)\p<a^FR ( 



(34) 



where Fj^ are as defined in Eq. and take again the 
integration over the perpendicular energy components 
into account. 

We solve the Wigner transport equation for the steady 
state with df /dt — numerically using the discretization 
scheme explained in Ref. 0- The position coordinate q 
is simply discretized as Nq equally-spaced points with 
the spacing Ag = Wn/{Nq — 1). Once A, is fixed the 
Fourier completeness relation gives the grid spacing Ap 



as a function of Aq as 



Pk 



TT 

a: 



k-i 



fc = l,2. 



(35) 



V{q,p) — 2 sm(pr) x 
Jq 

{Veff{q+lr)~Veff{q-^r)]dr.{31) 

is a one-dimensional potential kernel. The effect of Vef / 
is non-local, incorporating quantum interference effects. 
Namely, the potential kernel V{q,p — p') spreads accord- 
ing to Eq. (|30|l the WF f{q,p) among different values of 
p and adds the interference between alternative paths to 
the formulation. 

The electron density n^{q) and current density J{q) 
are calculated as 



1 r 



f{q,p)pdp 



(32) 



The derivative of ,f{q,p) with respect to q in the trans- 
port equation H3U|) has been calculated using the second- 
order differencing scheme (SDS). I.e., 



df{q,p) 



dq 



± 



Sf{q,p) - 4/(g ± Aq,p) + f{q ± 2Aq,p) 
2A„ 



(36) 



where the different signs are chosen according to the signs 
of ph in such a way that the proper boundary conditions 
given by Eq. 1)341) are coupled to the transport equation. 
It has been shown by Buot and Jensen 23, 24] that the 
SDS gives considerably better results than the conven- 
tional upwind-downwind differencing scheme (UDS). 
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D. Differences between the Green's function and 
Wigner function formalisms 

The GF and WF formalisms are related. Namely, we 
can first calculate the density operator in the GF formal- 
ism as 

/oc 
G"^ {x , x' ; uj) duj . (37) 
-oo 

Then the coordinate transformation of Eq. (|26() is per- 
formed and we get the WF f{p, q). If f{p, q) is calculated 
in this way Eqs. (|33|l and (|32|l of the WF formalism give 
exactly the same results as Eqs. ^ and (I19|l of the GF 
formalism. 

The differences in the results of the two formalisms are 
caused by the calculation of f{p, q) in the WF formalism 
by using the Liouville--von Neumann equation. For an 
infinite system the results would still be equal but a fi- 
nite calculation region causes the differences. Namely, 
the WF formalism does not have similar totally open 
boundary conditions as the GF formalism. Differences 
arise in this case also because there is no energy depen- 
dence in the Liouville-von Neumann equation whereas 
G< depends on oj. 

Using the GF formalism typically means that we have 
to calculate a lot of information which is of no further 
use. The calculation of G^'{x,x';lj) using N basis func- 
tions means inverting a N x N matrix. This requires the 
solution of N linear equations each having N unknown 
variables. The coefficient matrix is sparse, including only 
overlapping terms of the basis functions. In order to 
calculate the electron density we have to integrate over 
uj, which results in calculating G^{x,x';io) many times. 
Luckily, in the integrals n_2 and n_3 (Eqs. and 
(tTTjl l we need only two linear equations for G^{xl,x) 
and G^{xr, x). 

In the WF formalism we need to solve for the Wigner 
function, which depends on two variables p and q. Be- 
cause of the use of the discretization form of Eq. 
the dependence between these variables is more compli- 
cated than that between x and x' in the GF formalism. 
This means that we solve a set of linear equations which 
includes Np x Nq unknown variables. The coefficient ma- 
trix has a belt type filling where the width of the belt is 
large. This makes the linear equation hard to solve. 



IV. RESULTS 

In this section we give results of our electronic struc- 
ture and electron current calculations for the magnetic 
RTD structure described in Sec. IIA. We concentrate on 
the feasibility of the GF and WF formalisms to model 
this kind of systems. In the first subsection this is done 
by comparing the results of the two formalisms with each 
other. In the second subsection we compare our results 
with those measured by Slobodskyy et al. |l| and thereby 



analyze the electronic structure of the actual device. Fi- 
nally, we make predictions for the spin-polarized current 
when structure parameters of the magnetic RTD are var- 
ied. 



A. Comparison between the Green's function and 
Wigner function formalisms 

The electron density and the effective potential corre- 
sponding to the 0.15 V bias voltage and calculated using 
the GF and the WF formalisms are shown in Figs. 13 
and^ respectively. The two formalisms give very simi- 
lar results. There are Friedel oscillations in the electron 
density in the leads but their amplitude is so small that 
they are not visible on the scale used. The electron den- 
sity drops close to zero in the un-doped regions. For this 
value of VsD some of the low-energy resonance states 
are occupied. However, these states cause only a small 
density increment between the potential barriers. The 
electron density shows in both formalisms a small asym- 
metric behavior which becomes clearer with increasing 

VsD- 




20 40 60 80 100 120 

Position (nm) 



FIG. 3: Electron density calculated using the GF (solid line) 
and the WF (broken line) formalisms. The values of the tem- 
perature, magnetic field, and bias voltage used are 4.2 K, 6 T 
and 0.15 V, respectively. The gray areas denote level of posi- 
tive back ground charge. 

The effective potential in Fig. prises strongly at the 
interfaces between the doped and un-doped materials. 
Within the un-doped region the potential changes then 
rather linearly due to the applied bias voltage. The spin- 
splitting of the potential due to the external term of Eq. 
(0 does not propagate out of the quantum well region. 
Our GF and WF calculations do not include inelastic 
scattering. That would be important for finding self- 
consistent solutions if, on the higher potential side of 
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X (nm) Voltage (V) 



FIG. 4: Effective potential calculated using the GF (solid 
Une) and the WF (dashed line) formalisms. The values of 
the temperature, magnetic field, and bias voltage used are 
4.2 K, 6 T and 0.15 V, respectively. The effective potentials 
for spin-up and spin-down electrons differ remarkably only 
between the potential barriers. 



the barriers, there are "notch" states below the band of 
occupied electron states 0. In our case "notch" states 
do not appear because of the strong potential rise due to 
the un-doped material layers. 

The current through our magnetic RTD structure is 
shown in Fig. [Sjas function of the bias voltage Vsd- The 
spin-up and spin-down contributions, split strongly by 
the magnetic field of 6 T, show the typical RTD behavior. 
The most prominent peaks at around 0.14 V . . . 0.18 V 
actually correspond to the second lowest resonance states 
in energy. Most of the conduction takes place close to the 
right Fermi-level fj,ji because the large un-doped region 
diminishes strongly the tunneling probability at lower en- 
ergies. For this reason also the current peaks due to the 
lowest-energy resonances at around 0.03 V . . . 0.07 V are 
hardly visible in the GF results. They are clearer in the 
WF results, but this is partly due to numerical difficulties 
as will explain below. In the WF formalism the current 
may even change its direction and have negative values 
as seen in the lower panel of Fig. [5| 

The two resonance states per spin in the voltage region 
of . . . 0.25 V can be clearly seen in the local density of 
states (LDOS) calculated in the GF formalism. Fig. El 
shows the LDOS for the zero-magnetic field and the zero- 
bias voltage case. In the absence of the magnetic field, 
there is no difference between the spin-up and spin-down 
electron states. The first and the second resonance cor- 
respond to the quantum well states with zero and one 
node plane parallel to the layer structure, respectively. 
We can also see that especially at low energies the LDOS 
enhances just outside the un-doped region. 

The GF and the WF formalisms give almost equal elec- 



FIG. 5: Current through the magnetic RTD structure as 
a function of the bias voltage. The contributions due to 
spin-down (solid line) and spin-up (dashed line) electrons are 
shown. The results in the upper and lower panels are calcu- 
lated using the GF and the WF formalisms, respectively. The 
values of the temperature and magnetic field used are 4.2 K 
and 6 T, respectively. 




FIG. 6: LDOS of the GF formalism as a function of the 
position X along direction of the current and the energy u. 
The values of the temperature and magnetic field are 4.2 K 
and T, respectively. 



tron densities and effective potentials as demonstrated in 
Figs. 13 and 01 However, the requirements to reach simi- 
lar numerical convergences are very different. In the GF 
calculations we have used 36 fifth-order elements corre- 
sponding to 181 basis functions. A further increase of 
the number of elements does not change the results. In 
the WF calculations we have used as many as Nq = 320 
and Np — 300 discretization points. This means that the 



9 



number of unknown variables in the linear equations to 
be solved is of the order of 10^. If we reduce the number 
of discretization points from this magnitude, the effec- 
tive potential in the RTD region moves down becoming 
saggy. The very different numerical characteristics of our 
WF and the GF calculations demonstrate the efficiency 
of the high-order polynomial basis used in our GF imple- 
mentation. 

In the current- voltage curves of Fig. [S] the differences 
between the WF and GF formalisms are more visible 
than in the electron density of Fig. |3| although the for- 
malisms give the same qualitative behavior. The reso- 
nance peaks of the WF results are located at slightly 
lower bias voltages (energies) than those in the GF re- 
sults. This is in accordance with the lower effective po- 
tential in the WF calculations (See Fig. Ej). The GF 
formalism gives always by definition a positive current 
whereas the current calculated by using the WF formal- 
ism may become negative. This is a well-known artifact 
of the WF formalism. If we reduce the number of discre- 
tation points the current attains even more negative val- 
ues. This implies that the WF results of Fig. Oare not 
numerically fully convergent. However, with the present 
computer memory limits we cannot increase the number 
of discretization points much beyond 10^. At this stage 
the WF calculation begins to require also more CPU time 
than the GF calculations. This is somewhat surprising 
because one would think that the less approximative GF 
formalism would be computationally heavier. The use of 
a more accurate differencing scheme could improve the 
WF results without increasing the number of the dis- 
cretization points, but this would increase the filling of 
the coefficient matrix to be inverted and the CPU time 
needed. 

The numerical problems in the WF calculations are 
caused by the narrow resonance peaks. A Fourier trans- 
formation has to be evaluated over electron energies, 
and the resonance peaks require a remarkable increase 
of discretation points. Indeed, Fig. |31 shows that the 
current from the WF calculations for both spins oscil- 
lates strongly in the region of the first resonance peak, 
which is very narrow. The problem is not faced in the 
electron density calculations, because the contributions 
of the resonance states to the total electron density are 
small. Our real-space GF implementation can handle the 
narrow resonance peaks better because we use the adap- 
tive Simpson integral routine to calculate the density in- 
tegrals n_2, n-3, and the current integral along the real 
Lv axis and because the integral n_i is performed in the 
complex plane where the resonance peaks are broadened. 



B. Comparison to the experiments 

We compare our calculated results with the recent 
measurements by Slobodskyy et al. P|. In particular, 
we have calculated the current vs. bias voltage curves 
using the same magnetic field and temperature values as 



they have used. 




0.1 0.15 0.2 

Voltage (V) 




0.1 0.15 
Voltage (V) 



FIG. 7: Current through the magnetic RTD structure as a 
function of the bias voltage and different values of the mag- 
netic field. The structure is defined in Fig. Qand in Table 
I. The results in the panels a) and b) are calculated using 
the WF and the GF formalisms, respectively. The tempera- 
ture is 4.2K. For clarity, the successive curves are shifted by 
0.25 mA//im^ with respect to each other. 

The results of the WF and GF calculations for the cur- 
rent as a function of the bias voltage are shown in Fig.d 
corresponding to different magnetic fields and in Fig. |S1 
to different temperatures. As was mentioned above we 
have chosen the height of the potential barriers so that 
the widths of the two prominent spin-up and spin-down 
peaks from our GF calculations agree at low tempera- 
tures and high magnetic fields with experiment as well as 
possible. With this fitting only, our predictions for the 
magnetic field and temperature dependences are rather 
similar to experiments. The spin-up and spin-down con- 
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FIG. 8: Current through the magnetic RTD structure as a 
function of the bias vohage in different temperatures. The 
structure is defined in Fig. 0and in Table I. The results in 
the panels a) and b) are calculated using the WF and the GF 
formalisms, respectively. The magnetic field is 6T. For clarity, 
the successive curves are shifted by 0.25 mA//im^ with respect 
to each other. 



tributions separate in similar ways with increasing mag- 
netic field and merge together with an increasing tem- 
perature. These overall features are related to the de- 
pendence of the spin splitting Ai? on magnetic field and 
temperature in Eq. One should note that the dis- 

tances between the spin-up and spin-down peaks in the 
current are by a factor of two larger than AE (See the 
values in Table II). The reason is that part of the voltage 
is dropped over the first potential barrier. This fact was 
taken into account by Slobodskyy et al. by a lever 
arm of 2.1 when comparing the measured voltage split- 
ting with the theoretical AE. Our simulations confirm 



the magnitude of the lever arm. 

A more detailed comparison reveals that in the theo- 
retical current-voltage curves the resonance peaks are at 
about 30 mV higher voltages than in the experimental 
results. Several effects may contribute to this difference. 
First of all, uncertainty in the thickness of the quantum 
well affects strongly the resonance positions. The value 
of the series contact resistance in the measurements, the 
lack of inelastic scattering and the use of a constant ef- 
fective electron mass in our modeling may also shift the 
peaks. However, the energy shifts due to these uncer- 
tainties are expected to be smaller than the distances 
between the resonances for each spin direction: the dis- 
tance between the first and the second resonance is about 
100 mV (see Fig. (S)). Therefore our conclusion is that the 
current peaks seen in the measurements by Slobodskyy 
et al. correspond to the second resonance state in energy. 



C. Effects of the width of the potential well 

In the previous sections we have concentrated on mod- 
eling and analyzing the magnetic RTD structure by Slo- 
bodskyyei al. Next we predict how the geometry of the 
device affects the spin polarization of the current. The 
obvious parameter to be varied is the width of the quan- 
tum well because it determines the positions of the res- 
onances on the energy axis and affects strongly even the 
qualitative features of the current-voltage curves. Thus, 
in the following we keep the width and height of the po- 
tential barriers the same as in the previous calculations 
and vary the width Lss of the potential well using the 
values Ls5 = O.25L55, O.SLgg, and SL^^, where Lg^ is 
the original width of 9 nm (Table I). 

The current as a function of the bias voltage for the dif- 
ferent widths of the potential well is shown in Fig. O We 
see that the narrowing of the well moves resonance states 
towards higher voltages and the first resonance peak for 
each spin becomes active, i.e. its contribution to the to- 
tal current becomes evident. We also note that the mag- 
nitude of the current increases rapidly with decreasing 
potential well width. The increase of the potential well 
width to 2Ls5 causes the resonances to come closer to 
each other and overlap more strongly. This might be an 
undesirable feature for device applications. 

Another view of the functioning of the magnetic 
RTD structure as a spin switch is given in Fig. [TUI 
showing the relative polarization of the current, i.e. the 
difference between the spin-up and spin-down electron 
currents divided by the total current. For structures 
with very narrow potential wells the polarization is 
reduced, although the magnitude of the current is large 
(Fig. O. For structures with wide potential wells 
the spin splitting for a given magnetic field strength 
may become comparable with the distances between 
the resonances. This causes a complicated polarization 
structure as a function of the bias voltage as can be 
seen in Fig. ^| for the widest potential well of 2Ls5. In 
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FIG. 9: Current through a magnetic RTD structure as a func- 
tion of the bias voltage. The structure is defined in Fig. 
and in Table I with the exception that the width of the quan- 
tum well (region S5) is scaled relative to the original width 
I/gj. The results are calculated using the GF formalism and 
they correspond to the the magnetic field of 6 T and zero 
temperature. For clarity, the successive curves are shifted by 
0.35 mA/^m^. The lowest curve for 0.25xl/5g is scaled by a 
factor of 1/10. 
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FIG. 10: Relative polarization of the current through a mag- 
netic RTD structure as function of the bias voltage. The 
curves correspond to the total currents shown in Fig. |U] See 
also the caption of Fig. |0| 



conclusion, for a given magnetic field strength there is 
an optimum potential well width with good polarization 
properties and high current intensity. In the structures 
studied above and in the strong field of 6 T it would 
be approximately 0.3 . . . 0.5 times the original width ig^. 

V. CONCLUSIONS 

We have used the semiclassical Wigner function for- 
malism and the quantum-mechanical Green's function 
formalism to model spin-dependent electron transport 
through a magnetic resonant-tunneling diode structure. 
Our calculations are based on the self-consistent solu- 
tion of the electronic structures and currents within the 
density-functional theory. Both formalisms give very 
similar results for the electron density, but the current is 
more sensitive to the formalism used. We have traced the 
differences to be largely due to the numerical accuracy 
which is costly to achieve in the conventional implemen- 
tation of the Wigner function formalism. Surprisingly, 
the Green's function formalism implemented by using the 
finite-element method requires less computer resources 
and numerically converged results are obtained. 

Using the two schemes, especially the Green's function 
formalism, we have analyzed the recent measurements 
by Slobodskyy et al. |lj] for an actual magnetic resonant- 
tunneling diode. The magnetic-field and temperature de- 
pendencies in the measured current-bias-voltage curves 
are well reproduced. The two-peak structure is found to 
result from the spin-split resonance second lowest in en- 
ergy. We show that for a given strength of the magnetic 
field causing the spin-splitting there exists an optimum 
width for the quantum well so that high spin polarization 
and current intensity are achieved. 
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